Invention Journal of Research Technology in Engineering & Management (IJRTEM) 

ISSN: 2455-3689 

www.ijrtem.com Volume 3 Issue 1 1 January 2019 II PP 23-30 

An A-Stable Uniform Order Six Linear Multi-Step Method for Testing 
Stiffly First Order Ordinary Differential Equations. 

^abo John, 2 Domnic Raymond, 3 Adamu A. Bambur, 2 Yusuf T. Kyagya 

'Department of Mathematic, Adamawa State University, Mubi-Nigeria. 

2 Department of mathematics/Statistics, Taraba State Polytechnic, Suntai, Jalingo Campus-Nigeria. 
Department of mathematics/Statistics, Federal University Wukari, Taraba State. 

ABSTRACT: This paper formulates an A-Stable uniform order six linear multi-step method for testing stiffly 
first order ordinary differential equations. Power series approximate solution was used as the basic function using 
interpolation and collocation techniques to derive the continuous schemes which was found to converge. The 
developed method is then applied to solve the system of first-order stiffly ordinary differential equations and the 
result is displace side by side in tables with the existing once, the method tend to give better approximation when 
implemented numerically than the existing method with which we compared our results. 
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I. INTRODUCTION 

In the early 1950s, as a result of some pioneering work by (Curtiss, C. F. and Hirschfelder, J. O., 1952), (Stuart, 
A. M. and Humphries, A. R.,1996) was realized that there was an important class of ordinary differential equations 
(ODEs), which have become known as stiff equations, which presented a severe challenge to numerical methods 
that existed at that time. Since then an enormous amount of effort has gone into the analysis of stiff problems and, 
as a result, a great many numerical methods have been proposed for their solution. More recently, however, there 
have been some strong indications that the theory which underpins stiff computation is now quite well understood, 
and, in particular, the excellent text of (Hairer, E. and Wanner, G., 1996) has helped put this theory on a firm 
basis. As a result of this, some powerful codes have now been developed and these can solve quite difficult 
problems in a routine and reliable way. 

Interestingly, differential equations arising from the modeling of physical phenomena often do not have exact 
solutions. Hence, the development of numerical methods to obtain approximate solutions becomes necessary. To 
that extent, several numerical methods such as finite difference methods, finite element methods and finite volume 
methods, among others, have been developed based on the nature and type of the differential equation to be solved. 

A differential equation can be classified into Ordinary Differential Equation (ODE), Partial Differential Equation 
(PDE), Stochastic Differential Equation (SDE), Impulsive Differential Equation (IDE), Delay Differential 
Equation (DDE), etc. (Stuart, A. M. and Humphries, A. R.,1996). The main purpose of this paper is to outline 
some of the important theory behind stiff computation and to direct users of numerical software to those codes 
which are most likely to be effective for their particular problem. In sciences and engineering, mathematical 
models are formulated to aid in the understanding of physical phenomena. The formulated model often yields an 
equation that contains the derivatives of an unknown function. Such an equation is referred to as Differential 
equation, (Fatunla, S. 0..1988). As opposed to one-step methods, which only utilize one previous value of the 
numerical solution to approximate the subsequent value, multistep methods approximate numerical values of the 
solution by referring to more than one previous value (Skwame, Y., Sabo, J. and Kyagya, T. Y., 2017). 
Accordingly, multistep methods may often achieve greater accuracy than one-step methods that use the same 
number of function evaluations, since they utilize more information about the known portion of the solution than 
one-step methods do. The methods of Euler, Heun, Taylor and Runge-Kutta are called single-step methods 
because they use only the information from one previous point to compute the successive point, that is, only the 

initial point (/q , ) ! q ) is used to compute (t 1 , y 1 ) and in general Y[, is needed to compute .In this article, 
the system of first-order ODEs of the form 

y(x) = f(x,y(x)), y(x 0 )=y 0 (l) 

Is considered. This research is organized as follows: in the coming section, we carried out the derivation of the 
method, where we considered one-step through interpolation and collocation at first derivatives. 
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The details of the analysis of the method which include order, error constant, consistency, zero stability and 
stability region were discussed in Section three. In the fourth section, some numerical problems were solved and 
the performance of the developed method was compared with those of the existing methods, (Skwame, Y., 
Kumleng, G. M. and Bakari, I. A., 2017, Skwame, Y., 2018, Skwame, Y., Sunday, J. and Sabo, J. 2018). Finally, 
the conclusion was drawn in section five. 


II. THEORETICAL PROCEDURE 

In this section, the formation of single-step continuous uniform order six linear multi-step method for testing 
stiffly first order ordinary differential equations (1) is developed following (Alkasassbeh, M. and Zurni, O., 2017, 
Omar, Z. and Sulaiman, M., 2004). 

The general solution for solving (1) is derived in this section. Let the power series of the form 


v+m -1 


}’M= Ya U i 


i=0 


x-x n 

h 


( 2 ) 


be the approximate solution to (1) for X £ , X n+I ] where n = 0,1, 2, • ■ • N — 1, a's are the real 

coefficients to be determined, V is the number of collocation points, 1TI is the number of interpolation points and 
h = x„-x„_ i is a constant step size of the partition of interval [< 2 , b\ which is given by 
a = X Q < Jj < • • • < x N = b 

Differentiating (2) once yields: 
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Interpolating (2) at the point X n and collocating (3) at all points in the selected interval, i.e.. 


X n’ X 


1 , X 2 i x 31 X 4 and x n+l , gives the following nonlinear equations which can be written in 
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Applying the Gaussian elimination method on (4) gives the coefficient <2, S, foY 1—0 (l )l . 
These values are then substituted into (2) to give the implicit continuous hybrid method of the form: 
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y W= Yu a i Wy B+i +W/„ + i.+^ZA W/ b+ * 

!=0 1=0 i=k 


12 3 4 
_ 5’5’5’5’ 


On evaluating (5) yield the following continuous schemes, 

a 0 = 0 

B , = —f 2 /z(l44-616f + 1065f 2 -840f 3 +250 4 ) 

J 288 V 7 

/?, = — f z /t(72-428f + 885r -780r 3 + 250f 4 ) 

2 144 V 7 

/? 3 = - — r/t(-48-312/ + 735f 2 + 470f 3 ) 

^ 3 144 V ’ 

J3 4 = - — r/t(36-244/ + 615r -660 3 +250 4 ) 
- 288 

/?, = — r/t(50r + 4 Of + 9)(- 4 + 5f) 2 
288 


evaluate (6) at both grid and off-grid points yield the five continuous schemes as follows 
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(7) 

HI. CONVERGENCE ANALYSIS 

In this section, the basic properties of (7) shall be analyzed. 

Order and error Constants of the Methods 
Given linear difference operator (Butcher, J. C., 2009) 

t[y(x),h] = Y J [(d j y( x +jh)~ h J3jy\ x +jh)\ 

7=0 

and using Taylor expansion about the point x, we get 

£ [ y(x),h] = c 0 y(x ) + cjty 1 (*) + ••■ + ch q y q (x) + ■ ■ ■ 

where 


( 5 ) 


( 6 ) 


«•> 
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c 0 = a Q + a x + ••• + cc k 

c, = a x + 2 a 2 + —i- ka k — (/? 0 + /?, • • • J3 k ) 


c, = i (a, + Ta 2 + ■ ■ ■ + Va ,) - ^ (A, + 2'"" A + • • • + *" H 'A), <? = 2,3, • 

Such that when C Q = Cj = • ■ * = C p =0 and C p+ \ ^ 0 of (6) then 


r 

p+\ is the Error constant, and p is the order of LMM. 

According to (Butcher, J. C., 2009), the order of the new method is obtained by using the Taylor series and it is 
found that the developed method (7) is of uniformly order Six, with an error constants vector of: 


C 6 = [-1.8265x10 7 , -1.2529x10 7 , -1.6574xl0~ 7 , -1.0836x10 7 , -2.9101x10 7 f 
Consistency 

Definition 3.1: The hybrid block method (7) is said to be consistent if it has an order more than or equal to one 
i.e. P > 1. Therefore, the new method is consistent (Dahlquist, G., 1956). 


Zero Stability 

Definition 3.2: The hybrid block method (7) is said to be zero stable if the first characteristic polynomial ir(r) 

having roots such that |f, | <1 and if |r_| =1, then the multiplicity of T_ must not greater than two 
(Dahlquist, G., 1956 and Butcher, J. C., 2009). 

In order to find the zero-stability of hybrid block method (7), we only consider the first characteristic polynomial 
of the method according to definition (3.2) as follows 
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which implies r = 0, 0, 0, 0, 1. Hence the method is zero-stable since G | < 1 


and if 



= 1 . 


Convergence 

Theorem (3.1): Consistency and zero stability are sufficient condition for linear multistep method to be 
convergent. Since the method is consistent and zero stable, it implies the method is convergent for all point 
(Butcher, J. C., 2009). 

Regions of Absolute Stability (RAS) 

Definition 3.3: The Region of Absolute Stability is the set of all points Z £ Z such that all roots of characteristic 
equation are of absolute value less than one, (Lambert, (1973). According to (Dahlquist, G., 1956 and Lambert, 
J. D., 1991) the absolute stability region of the new method is shown in the figure below. 
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F=ie(z) 

Figure 3.1: Region of Absolute Stability 

IV. THE IMPLEMENTATION OF METHOD 

In this section, the efficiency and the performance of the new method (7) is investigated with four systems first 
order stiffly ODEs. The problems considered are the ones widely used by (Skwame, Y., Sabo, J. and Kyagya, T. 
Y., 2017, Sabo, J., Skwame, Y., Kyagya, T. Y., (2018), Akinfenwa, O. A., Abdulganiy, R. I., Akinnukawe, B. I., 
Okunuga, S. A. and Rufai, U. O., 2017), and the results obtained shall be compared with the existing method in 
terms of error. 


Problem Consider the stiffly problem, 
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2e 4111 (sin(4(k)- cos(40x)) J 

(Source: Skwame, Y., Sabo, J. and Kyagya, T. Y., 2017) 


*4 


O < jc > 0.8, 


h = 0.1 


Table 4.1: Comparison of result of new method with that of (Skwame, Y., Sabo, J. and Kyagya, T. Y., 

2017) 


x-values 

Error in Skwame, et-al., (2017) 

Error in new method 

Y, 

y 2 

^3 

Yi 

a 2 

y 3 

0.1 

2.23xl0" 2 

2.23xl0" 2 

2.53xl0" 2 

2.15xl0" 2 

2.15xl0" 2 

2.42x10 -2 

0.2 

1.06xl0” 4 

9.14xl0" 5 

1.68xl0" 4 

9.65 xlO" 5 

80 

0*8 

X 

o 

1 

80 

bo 

X 

o 

0.3 

8.23xl0” 6 

9.10xl0“ 6 

1.33xl0” 5 

4.22x10” 7 

4.21xl0” 7 

1.24xl0” 7 

0.4 

9.60xl0" 6 

9.30xl0" 6 

1.60xl0” 7 

1.30x10 7 

1.30x10 7 

1.35xl0 " 7 

0.5 

9.67x10 6 

9.67x10 6 

cp 

o 

X 

oo 

so 

6 . 00 x 10 10 

2 . 00 x 10 10 

5.89x10 10 

0.6 

9.50x10 6 

9.50x10 6 

9.12x10 11 

4.00x10 10 

0 . 00 x 10 ° 

3.75x10 11 

0.7 

9.08x10 6 

9.08x10 6 

1.05x10 10 

4.00x10 10 

0 . 00 x 10 ° 

2.59x10 10 

0.8 

8.49x10 6 

8.49x10 6 

6.84x10 11 

4.00x10 10 

0.00x10° 

1.30x10 10 
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Problem 4.2 

Consider the stiffly problem, W = —lOOv! +9.901 y 2 ; y,(o)= 1 

y 2 = 0.1 y l -y 2 ; y 2 (o) = 10, h = 0.1 

With Exact Solution 
.V, CO = 
y 2 (x)=]0e 099x 

x e [0, l] 

(Source, Sabo, J., et-al. 2018) 

Table 4.2: Comparison of absolute error with that of Sabo, et-al., (2018) 


x - values 

Error in Sabo, et-al., (2018) 

Error in new method 

>’i (*) 

y 2 (x) 

wC*) 

y 2 (x) 

0.1 

1.80xl0~ 9 

1.40x10 8 

5.00x10 10 

4.00 xl0~ 9 

0.2 

2.70x10 9 

2.30x10 8 

1.50xl0’ 9 

7.00 xlO^ 9 

0.3 

3.70x10 9 

3.30x10 8 

1.50x10 9 

9.00x10 9 

0.4 

4.40 xl0“ 9 

3.90x10 8 

2.20xl0~ 9 

1.30x 10~ 8 

0.5 

5.00x 10 9 

470x1O 8 

2.20xl0" 9 

1.50x10 8 

0.6 

5.20xl0" 9 

5.00x10 8 

2.30x 10~ 9 

1.30x10 8 

0.7 

5.40x 10” 9 

5.20x 10~ 8 

1.30x10 9 

1.20x10 8 

0.8 

5.70 x 10“ 9 

5.40x10 8 

2.10xl0~ 8 

1.40x10 8 

0.9 

5.60x 10~ 9 

5.50xl0~ 8 

1.70x10 9 

1.50x10 8 

1.0 

5.70 x 10“ 9 

5.50x10 8 

2.40x1 O' 9 

1.70x10 8 


Problem 4.3 

Consider the stiffly problem, 

y'i = -yi> > 7 i(°) =1 

y 2 = -2000y 2 ; y 2 (o) = l 

y l (x) = e~ x 

y 2 (x)=e-™* 


0.1, 0<X<1 


Table 4.3: Comparison of absolute error 


x - values 

35 


0.1 

2.00x10 10 

7.96x10 1 

0.2 

1.00x10 10 

6.34x10 1 

0.3 

3.00x10 10 

5.04x10 1 

0.4 

5.00x10 10 

4.01x10 1 

0.5 

5.00xl0~ 10 

3-19x10 1 

0.6 

5.00x10 10 

2.54x10 1 

0.7 

5.00xl0 1 ° 

2.02x 10 1 

0.8 

5.00xi(r 10 

1.61x10 1 

0.9 

5.00x10 10 

1.28x10 1 

1.0 

4.00x10 10 

1.02x10 1 
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Problem 4.4 Consider the stiffly problem, 

y 1 +0.1y 1 +49.9y 2 =0; y 1 (o) = 2 

y' 2 + 50y 2 = 0; y 2 (0)=0, h = 


y 3 -70y 2 + 1 20 y 3 = 0 ; y 3 ( 0)=2 

With Exact Solution 

yi (x) = e- Xx +e X,x 
>> 2 (x) = c 

y,(x) = e- 50 ' + e- l20x 

X G [O, l] 

(Source, Akinfenwa, O. A., et-al. 2017) 


Table 4.4: Comparison of absolute error 


x - values 

Absolute error 

Jr 


y 3 

0.1 

9.85x10 3 

5.00x10 10 

8.87x10 1 

0.2 

1 . 96 xl 0" 2 

1 .00x10 9 

7.40x10 1 

0.3 

2.93x10 2 

1.40x10 9 

6.98x10 1 

0.4 

3.88x10 2 

1.70x10 9 

6.19x10 1 

0.5 

4.83x10 2 

2 .00x10 9 

5.49x10 1 

0.6 

5.76x10 2 

2.30x10 9 

4.87x10 1 

0.7 

6.70x10 2 

2.50x10 9 

4.32x10 1 

0.8 

6 . 61 x 10 2 

2.70x10 9 

3.83x10 1 

0.9 

8.52x10 2 

2.90x10 9 

3.40x10 1 

1.0 

9.42x10 2 

3.10x10 9 

3.01x10 1 


V. CONCLUSIONS 

The single-step continuous uniform order six linear multi-step method for testing stiffly first order ordinary 
differential equations has been developed in this paper. The analysis of the new method was studied and it was 
found to be consistent, convergent and zero-stable, with the region of absolute stable within which the method is 
stable. The newly constructed method was applied to solve four systems of first-order stiffly ordinary differential 
equations and from the results obtained, it was clear that the developed method performed better. Therefore, the 
general solution of first order Linear Multistep Method (LMM) is a convenient technique for determining the 
solutions of mathematical modeling since it can approximate the result even though the efficiency is less than the 
other multistep method. 

Recommendation: The pair of numerical methods developed in this paper is recommended for testing first order 
stiffly ordinary differential equations. The basis function (power series) used is also recommended for the 
derivation of Numerical methods for first order differential equations and the pair of methods derived are also 
recommended for the solution of systems of first order stiffly ordinary differential equations. 
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